Strong coupling theory for the superfluidity of Bose-Fermi mixtures 
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We develop a strong coupling theory for the superfluidity of fermion pairing phase in a Bose-Fermi 
mixture. Dynamical screening, self-energy renormalization, and pairing gap function are included 
self-consistently within the adiabatic limit (i.e. the phonon velocity is much smaller than the Fermi 
velocity). An analytical solution for the transition temperature (T c ) is derived within reasonable 
approximations. Using typical parameters of a 40 K- 87 Rb mixture, we find that the calculated T c 
is several times larger than that obtained in the weak coupling theory, and can be up to several 
percents of the Fermi temperature. 



Recently superfluidity of fermion pairing phase has be- 
come an extensiv ely s tudied subject in the field of ul- 
tracold atoms 12 3 4 The attractive interaction 

for two fermions to form a Cooper pair can be provided 
either by the direct s-wave scattering between different 
spin states |jH3- or by the effective interaction medi- 
ated by the condensate phonons in a Bose-Fermi mix- 
ture (bfm) m. 

The former systems attract a lot of 
attention on the BCS-BEC crossover near the Feshbach 
resonance regime 0|, while the latter systems, which 
can mimic some traditional electron-phonon systems in 
solid state physics, are expected to have many interesting 
many-body phases in different parameter regimes 0Q@ • 

To observe the predicted fermion pairing phase or other 
many-body states in a BFM, experimentalists usually 
have to increase the boson-fermion interaction strength 
via a heteronuclear Feshbach resonance 00- However, 
by doing so the system will be driven into a strong cou- 
pling regime, where the direct application of the Bardeen- 
Cooper-Schrieffer (BCS) theory for the superfluidity of 
fermion pairing phase becomes unjustified and not re- 
liable. For example, the strong fermion-boson interac- 
tion may renormalize the effective mass of fermions and 
hence change its density of states near the Fermi sur- 
face. Besides, the fermionic quasi-particle fluctuations at 
finite temperature can also dynamically screen the effec- 
tive interaction induced by condensate phonons. None 
of the these effects have been included in the weak cou- 
pling theory developed in the literature |4j|5||fJ]. More- 
over, in most of the weak coupling theories, one usu- 
ally applys an instantaneous approximation (i.e. assum- 
ing the phonon velocity, q,, is much larger than the 
Fermi velocity, Dp). Such approximation is, however, 
very unrealistic in typical systems like 40 K- 87 Rb or 6 Li- 
23 Na mixtures, which are usually in the adiabatic regime 
(cb/vF ~ 0.05 — 0.5 < 1) due to the weak interaction 
between bosonic atoms |5J. Therefore in order to have 
a reasonable comparison between the theoretical results 
and the experimental measurement, it is necessary and 
important to develop a full strong coupling theory, in- 
cluding the adiabaticity of condensate phonons as well as 
the strong correlation effects of a BFM self-consistently. 

In this Letter, we develop a strong coupling theory for 
the superfluidity of (s-wave) fermion pairing phase in a 
BFM within the adiabatic regime by generalizing the cel- 



ebrated Migdal-Eliashb erg equ ations used in the conven- 
tional superconductors [l(illri| . By strong coupling, we 
mean to treat the fermion self-energy, dynamical screen- 
ing of the effective interaction, and the pairing gap func- 
tion self-consistently. Hence the obtained results are re- 
liable even when the fermion-boson interaction is not 
weak. Vertex diagrams are safely neglected by apply- 
ing Mi gdal's theorem in the adiabatic regime of phonons 
|lfjllj . Using a single mode approximation we further 
derive an analytical solution for the superfluidity tran- 
sition temperatures and compare them with the known 
weak coupling BCS results in different parameter regime. 
In a typical 40 K- 87 Rb mixture, we find that the obtained 
T c is several times larger than the weak coupling results, 
and can be up to a few percents of the Fermi energy 
before reaching phase instability. Effects due to strong 
boson-fermion correlations are also critically discussed. 

We consider a three dimensional BFM composed of 
spin polarized bosons and fermions in two equally occu- 
pied hyperfine spin states. For simplicity, we neglect the 
inhomogeneous magnetic trap potential so that the total 
Hamiltonian can be written to be [a 



H = 



(1) 



^Plp b - k + Unfold + I nii i> J k .. 



where 6k and /k, s are the field operators for bosonic and 
fermionic atoms of momentum k and isospin s =t)|- 
kp+k^p * s the boson density operator, and — 
k 2 /2m{, — fib is the bosonic kinetic energy with nib being 
the atom mass and fib being the chemical potential. Simi- 
lar notations also apply to fermions with super/subscript 
/, and ^ ee pl A + p[ d . U bb = U„ = 



and 



are respectively the boson-boson, fermion- 



fcrmion and boson-fermion pseudo-potential strength, 
with dij being the associate s— wave scattering length. 
m r is the reduced mass and VI is the system volume. 

Throughout this Letter we assume that the tempera- 
ture is so low that all bosons are condensed at zero mo- 
mentum state. The condensate excitations (phonons) are 
of Bogoliubov type dispersion: w£ = Cfc|k|-y/l + |k| 2 £ 2 , 
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FIG. 1: Feynman diagrams studied in the strong coupling 
theory of a BFM. (a) Full (dressed) effective interaction (see 
the text); (b) Migdal-Eliashberg equations, where S and W 
are the self-energies for normal and anomalous Green's func- 
tions respectively; (c) Polarizability, P, is replaced by bubbles 
without vertex correction; (d) Vertex correction neglected in 
the strong coupling theory due to Migdal's theorem. 




FIG. 2: Phonon dispersions including (solid line) and not- 
including (dashed line) the dynamical screening effects. Be- 
low the dotted line is regime of quasi-particle(fermion) ex- 
citations, where the Landau damping broadens the phonon 
spectral function as shown in the inset. Inset: the real(solid 
line) and imaginary (dashed line) parts of the full effective in- 
teraction, Vjg(q,tj), at q — kp. Parameters are in Ref. [T^. 



where Cf, = s/nbUbb/mb and £b = \J '1 / 'AmbUbUbb are re- 
spectively the phonon velocity and the healing length 
of the condensate with density n^. As a result, one 
can effectively describe the BFM system by +bo fal- 
lowing fermion-phonon type Hamiltonian 4 5 6 



E k , s 4fLk + E k ^ k /3 k /3 k + ffk(/3 k + £! k )/£, 





: H = 




where 



/3k is the phonon operator and gu = UbfynbC^/cu^ 
measures the fermion-phonon coupling strength. In- 
tegrating out the phonon field, one can obtain a re- 
tarded phonon-induccd interaction between fermions: 

Vph(k,w) = ^p^a . When the phonon velocity 

is smaller than the Fermi velocity as considered in this 
Letter, the dynamical screening due to fermion density 
fluctuations has to be included by calculating the Dyson's 
equations shown in Fig. [Tfa), where U" — V p h and 
U = Vph + Uff are the bare effective interaction be- 
tween fermions in the spin parallel/perpendicular chan- 
nels. The full dressed effective interaction can then be 
easily derived to be 



V. 



ctr 



(1 _ u^XpyjhX + u ± PU A 
(l-U\\P) 2 - (PU^) 2 



(2) 



where -P(k, lu) is the polarizability due to fermionic quasi- 
particle and Cooper pair fluctuations. In Fig. |3we show 
a numerical result of V^(q,uj) and the screened phonon 
dispersion (obtained by tracking the peak position of 
lmV c jf(q,uj)) at T = T c for a typical 40 K- 87 Rb system 
fl2| . As we will see later, although the screened phonon 
dispersion is close to the bare (unscreened) results, their 
spectral weights can still be quiet different due to the 
fluctuations near the Fermi surface. 

Now the single particle Green's function, G(p, t) = 
— i\f/ p CT (*)/p CT (0)}, and the anomalous Green's func- 



tion F(p,t) = -i(f/_ Pji (i)/ P)T (0)> (where t is the 
time-ordering operator) can be calculated within a self- 
consistent Hartrec-Fock approximation as shown in Fig. 
H^b). Such meanfield treatment is justified in strong 
fermion-phonon interaction regime because the vertex 
correction is proportional to the ratio of the phonon ve- 
locity to the Fermi velocity (known as Migdal's theo- 
rem an d hence its leading order correction to the 
bare vertex (Fig. dd)) can be rather small (Fi < 0.1, 
see Ref. 12]) in typical 40 K- 87 Rb or 6 Li- 23 Na mix- 
tures. By the same reason, we can also approximate 
the full fermion polarizability, P, by a bubble diagram, 
P (k, W ) = |jEp/ %G(p, u)G(p + k, v + u), without 
vertex lines inside the bubble diagram (Fig. (He)). The 
contribution coming from pair fluctuations can be ne- 
glected since here we are only interested in the calcula- 
tion of T c . 

To formulate the Migdal-Eliashberg equations (MEE, 
see Fig. ^b)) for a BFM system, we note that the phonon 
excitation energy, is unbounded in the high energy 
limit, and hence the involved momentum integration can- 
not be restricted to the vincinity of the Fermi surface, 
as done in the original theory for the conventional super- 
conductors [llj . We therefore define the following energy 
dependent interaction strength: 



FS 



rf 2 p J B||,_ L (p - y/2m f Ez, v) 
v F (2n) 3 N(E F )/N(E) 



(3) 



where Bii j_(q, v) = 1 ImV l }\ 



weight of the dressed interaction; j FS 



(q, v) is the spectral 
denotes the in- 
tegration over Fermi surface, where the gap function is 
measured. N(E) = m f yj2mjE / '2tt 2 is th 3D density of 

states for E > 0. Note that a 2 F^ t x(v, E) drops to zero 
rapidly for E — > oo due to the sharp spectral function, 
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B\\ x, near the phonon excitation energy (see Inset of Fig. 
|2J). Following the same argument as Eliashberg ^l], we 
can neglect the momentum dependence of self-energies, 
S and W, and then solve Fig. ^b) by introducing the 
gap function, A(w) = W (u>) / Z (u>) , and a renormaliza- 
tion coefficient, Z(uj) = 1 - (S(uj) - S(-uj))/2uj [ll|. 



Z(ui — > 0) _1 < 1 measures the single particle spectral 
weight at Fermi surface, and is equal to one for nonin- 
teracting fermions. Details of the derivation for a BFM 
system will be presented elsewhere 13] . The obtained 
MEE are the following: 



LU{1 - Z{i0)) 

A(lu)Z(lu) 



dv 



de a 2 Fi,(>,/! + eZ(e))Re 



Uff J den F (e)Re 

-oo 

dv 



A(e)sgn(e) 



Ve 2 ~ A(e) 2 



^ 2 - A(e) 2 
dea 2 F±(v,n + eZ(e))Re 



+ eZ(e)) 



A(e)sgn(e) 
V? 2 ~ A(e)2 



n_F(e)+n s (^) 1 - nj? (e) + n B {v) 



— lu — v — iS 



e — u> + v — id 



n F (e)+n B (u) 1 - n F (e) + n B (v) 



— v — iS 



e — u) + v — % 



(4) 



(•5) 



where n F / B (e) — (e c / T ± is the Fermi/Bose 

distribution function. Z(e) is the real part of 
Z(e)y/l — (A(e)/e) 2 . Note that the instantaneous in- 
teraction, Uff, contributes to the self-energy, S, via a 
constant shift only [111 ] and therefore does not appear in 
the r.h.s. of Eq. (0}. Furthermore, the ultraviolet diver- 
gence caused by the pseudo-potential Uff has also been 
removed in Eq. J5J) by eliminating the bare Lippmann- 
Schwinger equations in vacuo as done in the weak cou- 
pling theory |2l4j . Therefore the MEE shown above have 
included both the non-retarded direct interaction and the 
retarded induced interaction on equal footing. 

To calculate the T c of Cooper pairs we first linearize the 
above equations by taking A(uj) — > and approximate 
Z{lS) by Zq = Z(u> — > 0). Assuming T c is much smaller 
than Fermi energy and typical phonon energy, we can 
simplify Eqs. ©-© to be 

i-*_- a r*r*^ii^M (6 ) 

Jo Jo + v) 2 

r°° de 

Z = Uff —n F (e)N(n + eZ ) 

J —oo e 



,/„ / Itfi^loflF^^ + Zoe). (7) 
-oo e(e - v) 



To solve T c and Zq 
ther apply the single 



analytically we fur- 
mode approximation, 



-7r ^im i ( i)fi( v ~ where 

11,-L/ 



^imC?) = ~ \ Io° IvoV^ (q, v) is the total spectral 
weight. The phonon energy, u> q , is obtained from the 

peak position of ImVdJf" L (g, v)- Within such approxima- 
tion a 2 FM_L of Eq. © can be greatly simplified and T c 
of Eq. Q can be solved to be 



E F 



8o 



7re 2 Zn 



exp 



l + A||(Z ) + Aio 8 (£>) 
— X±(Z n ) 



(8) 



where 7 ~ 1.781. A = N(E F )U ff and X±(Z ) = 

2 J^ ^a 2 F\(v, 0) are respectively the direct and in- 
duced coupling strengths for Cooper pairs. D = 
AE F /e 2 Zo is an energy scale of Fermi energy. Zq is then 
determined by solving Z n = 1 + X\\(Z ) self-consistently 
since Zq appears in the effective interaction via the po- 
larizability, P. The two coupling strengths, Aii(Zo) and 
Xiog(D), are given by (here v c = uj 2 k F ) 



\\{Zo) 



dvA^ (1 



v v + |e_(z/)| v + e+(y) 



dvA^ {v) 



v + e-{v) v + e+(u) 



'"^ D 2 {v+\e- \)(u + e+) 



dvA x (u) (e + {v + e-) 



log 



(9) 
(10) 



r 



where A" :J -(^) 



N(E F ) 
2kl 



duj q (q,) 
dq 



1 1*V-^{q*) with Q*( v ) g iven b y = v > and e ±(^) = ll/ 2m f ± 
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FIG. 3: (a) Superfluid temperature as a function of boson- 
fermion interaction strength, |af,/|fcF, for the strong cou- 
pling (solid lines) and weak coupling (dashed lines) theo- 
ries. Upper/lower pairs of lines are for boson density rib = 
Other parameters are given in Ref. PJ. (b) 
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Zq — 1 and — Aiog for the inset) in the same 



calculation of (a) with nt, = 100 um 



q* kp/rrif. 

It is instructive to compare the analytical T c in Eq. 
|jHJ with the weak coupling results 4J: one can see 
that Eq. (JHJ reproduces the weak coupling result Q: 



->0 _ 87 

7re 2 



Tx = -f^exp 

_ N(E F )UZ f 



An-A'l 



if |A||(Z )|,|Ai og (£>)| « 1. 



Here 



= // Voh'^'U m (l + (2fcF^f)) 2 ) is the phonon-induced 
coupling strength without any dynamical screening Q. 
(Note that here we only discuss the s-wave pairing sym- 
metry. The non s-wave scattering strength is in general 
very small and can be neglected.) In Fig. OJa) we com- 
pare numerical results of T c and T®, using typical param- 
eters of 40 K- 87 Rb systems [l2|. We find both of them 
increase sub-exponentially as the interaction strength, 
\abf\kp, increases. T c is larger than Xj? by a factor of 



2-8 for |a&/|fc,F < 0.12, and then becomes almost equal to 
Tj? for stronger interaction. Note that T c /T® is large for 
small \dbf\kp simply because the repulsive direct interac- 
tion, Uff > 0, sets the minimum value of \at,f \kp for pair- 
ing phase. This ratio will decrease to one as \abf\kp — > 
if Uff < 0. When \abf\kF > 0.15 the screened phonon 
velocity, q,, becomes imaginary and the system becomes 
unstable toward collapse. This critical value of \abf\kp 
is a little higher than the results predicted by the weak 
coupling theory Q due to the reduction of the density 
of states in the strong coupling theory. Analyzing each 
quantity of Eq. JSJl more carefully (see Fig. Efb)), we 
find that although the ratio of A^ to A^ increases as 
l^t/I^F increases via dynamical screening, the final ratio 
of T c to Xj? still becomes smaller due to the strong quasi- 
particle renormalization, Zq (see the inset), in the strong 
coupling regime. Numerically solving the full strong cou- 
pling equations (Eqs. Q-©) can also provide more in- 
formation about the gap function spectroscopy, A(w), 
and fermion single particle spectral function etc., which 
can be measured by using rf-spectroscopy and/or Bragg 
scattering spectroscopy in the present experimental tech- 
nique. 

In summary, we have derived the full strong coupling 
theory for the fermion s-wave pairing phase in a Bose- 
Fermi mixture. Our results apply to the limit of slow 
phonon velocity and hence are valid in the existing BFM 
system ( 40 k- 87 Rb or 6 Li- 23 Na). The predicted critical 
temperature for a typical 40 K- 87 Rb mixture can be as 
high as a few percents of the Fermi energy, which should 
be acheivable by present experimental groups. 

We thank fruitful discussion with J.H. Bao, E. Demler, 
M.D. Lukin, and C.-Y. Mou. 
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